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Four-dimensional Simulation of the Hot Electroweak Phase Transition 
with the SU(2) Gauge-Higgs Model * 
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'^Center for Computatfonal Physics, University of Tsukuba, Ibaraki 305, Japan 

We study the finite-temperature phase transition of the four- dimensional SU(2) gauge-Higgs model for inter- 
mediate values of the Higgs boson mass in the range 50<m£f <100GeV on a lattice with the temporal lattice size 
Nt = 2. The order of the transition is systematically examined using finite size scaling methods. Behavior of the 
interface tension and the latent heat for an increasing Higgs boson mass is also investigated. 



1. INTRODUCTION 

The possibility that the baryon number asym- 
metry of the Universe has been generated in the 
course of the electroweak phase transition has 
led to recent lattice investigations of the SU(2) 
gauge-Higgs model by several groups Their 
studies have shown that a first-order phase tran- 
sition takes place at a finite temperature in this 
system for a light Higgs boson, which, however, 
becomes rapidly weak as the Higgs boson mass 
increases. A central question, relevant for the 
baryon number asymmetry problem, is whether 
the first-order transition survives with a sufficient 
strength for realistically large Higgs boson mass, 
experimentally bounded by Mh > 64GeV[H. 
To answer this question many studies have been 
done within the perturbatively reduced three- 
dimensional model0]. On the other hand, only 
a few studies with the original four-dimensional 
model exist in this region of Higgs boson mass[|l|. 
In this article, we report results of our simula- 
tion of the four-dimensional model, aiming at a 
systematic finite-size scaling analysis of the order 
of the transition for the Higgs boson mass in the 
range 50<mff <100GeV. 

2. SIMULATION 

We employ the standard action given by 

^ = EiEfTrC/„,^. + 5]«;Tr($tC/„,M<i>„+A) 
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with the complex 2x2 matrix <i> decomposed 
as = PnOtn, Pn > 0; a„ G SU(2). All of 
our simulations are made for the temporal ex- 
tent Nt = 2. We set the gauge coupling (3 = 
— 8, and make simulations for 6 values of 
the scalar self-coupling A as listed in Table [|. 
Also listed in the table are the zero-temperature 
Higgs boson mass Mh at the transition point of 
Nt — 2 lattice, estimated by interpolating avail- 
able data for M//|^-||] as a function of A. For 
each value of A runs are made on an N'^ x 2 lattice 
with Ns = 8, 12, 16, 24, 32, and in addition with 
N, = 40 for A = 0.001, 0.0017235, 0.003(Af// = 
67, 85, 102GeV ). Gauge and scalar fields are up- 
dated with a combination of the heat bath|^] and 
overrelaxationl^ algorithms in the ratio reported 
to be the fastest in ref. H). For each parameter 
point we make 10^ iterations of the combined up- 
dates. 

3. FINITE-SIZE SCALING ANALYSIS 

Let us define the angular part of the spatial 
component of the hopping term in the action by 
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In Fig. |l| we show the volume dependence of the 
maximum height of the susceptibility defined by 



(3) 
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with V ~ N^, which is calculated by the stan- 
dard reweighting technique. For 47 < Mh < 



2 



Table 1 



Choice of scalar self coupling and corresponding zero temperature Higgs boson mass. 



A 


0.0005 


0.000625 


0.00075 


0.001 


0.0017235 


0.003 


Mjj(GeV) 


47 


53 


58 


67 


85 


102 




60000 



Figure 1. Maximum height of the susceptibility of 
As as a function of the volume. Lines are guides 
for eyes. 
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Figure 2. Valley depth of Binder cumulant of 
As as a function of the inverse volume. Filled 
symbols show the infinite volume limit. 



67GeV the maximum value increases linearly to- 
ward larger volumes, which is expected for the 
case of a first-order transition. In contrast, we 
observe a very flat volume dependence for 85 < 
Alfj < 102GeV, albeit the maximum value is in- 
creasing slowly in the range of volume used here. 

In Fig. H we plot the valley depth of the Binder 
cumulant defined by 



1 _ iS^'s) 



(4) 



3(A2)2- 

as a function of inverse volume. Lines are lin- 
ear fits to the largest 3 volumes for each Mh- 
For 47 < Mh < 67GeV the value extrapolated to 
the infinite volume clearly deviates from 2/3, pro- 
viding additional evidence for a first-order tran- 
sition. For 85 < AIh < 102GeV the deviation de- 
creases by an order of magnitude, although still 
finite within the error. 

These results clearly show that the transition 
is of first order for 47 < Mh < 67GeV. It is 
also clear that the transition, if first order, is a 
very weak one at Mh = 85GeV and 102GeV. It 
is quite possible that the transition has turned 



into a crossover for this range of Mh- Data for 
larger volumes are needed, however, for a conclu- 
sive analysis on this point. 

4. LATENT HEAT AND INTERFACE 
TENSION 

The Higgs boson mass dependence of the la- 
tent heat Ae provides a physical measure of the 
weakening of the first-order transition as Mh 
increases. Here we calculate Ae for Mh = 
47, 53, 58, 67GeV using the Clausius-Clapeyron 
equation M, 



Ae ~ -MfjKA(p^ 



(5) 



where the gap A(p^) is estimated from the run 
on the largest volume for each Mh- The result is 
shown in Fig. |[ 

Interface tension provides another indicator of 
the strength of the first-order transition. Let 
Pmax and Pmin bc the peak and valley height 
of the distribution of As, reweighted such that 
the two peaks have an equal height. Define 
av = -{N^/2N^)\ii{Prmn/P,nax)- For Spatially 
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Figure 3. Latent heat as a function of Higgs boson 
mass. Open symbol is from ref. [4] and is shifted 
shghtly left for visualization. 

cubic lattices, finite-size formula for the true in- 
terface tension a is given by|^ 
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where c is a constant independent of Ng . Making 
a two parameter fit of ay obtained for the largest 
three volumes for Mh — 47, 53, 58GeV or two 
volumes for Mh — 67GeV we find cr/T^ shown in 
Fig. I 

Both the latent heat and the interface tension 
rapidly decrease with an increasing Higgs boson 
mass and seem to vanish around Mh ~ 80GeV. 

5. SUMMARY 

Our finite-size scaling study establishes a first- 
order transition for 47 < Mh < 67GeV. For 
larger Higgs boson masses a rapid weakening of 
the transition makes it difficult to draw a defini- 
tive conclusion on the order within the range of 
lattice volumes employed in our simulation. How- 
ever, combining finite-size data with results for 
the latent heat and the interface tension, our four- 
dimensional study suggests that the first-order 
transition terminates around Mh ^ 80GeV in 
the Nt = 2 SU(2) gauge-Higgs model. This is 
consistent with the result of a recent finite-size 
scaling study carried out in the dimensionally re- 
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Figure 4. Interface tension as a function of Higgs 
boson mass. Legends are the same as in Fig. 3. 

duced three-dimensional model pO|. 

ACKNOWLEDGEMENTS 

I would like to thank Akira Ukawa for useful 
discussions. The numerical calculations were car- 
ried out on VPP/500 at Science Information Pro- 
cessing Center of University of Tsukuba and at 
Center for Computational Physics at University 
of Tsukuba. 

REFERENCES 

1. For reviews see K. Jansen, Nucl. Phys. B 
(Proc. Suppl.) 47 (1996) 196; K. Rum- 
mukainen, in these proceedings. 

2. P. Janot, Nucl. Phys. B (Proc. Suppl.) 38 
(1995) 264. 

3. B. Bunk et al, Phys. Lett. B 284 (1992) 371. 

4. Z. Fodor et al., Nucl. Phys. B 439 (1995) 147. 

5. F. Csikor at al, |hcp-lat/9601016i 

6. B. Bunk, Nucl. Phys. B (Proc. Suppl.) 42 
(1995) 566. 

7. Z. Fodor and K. Jansen, Phys. Lett. B 331 
(1994) 119. 

8. W. Buchmiillcr, Z. Fodor and A. Hebecker, 
Nucl. Phys. B 447 (1995) 317. 

9. Y. Iwasaki at al, Phys. Rev. D 49 (1994) 
3540. 

10. K. Kajantic at al, |hep-ph/9605288 . 



